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ABSTRACT 

Theoretical and experimental investigations were conducted for the 
purpose of improving accuracy in calculating pressure distributions for 
various tension shell decelerator shapes and evaluating the effects on 
derived shapes resulting from the application of different pressure dis- 
tributions. Although tension shell shapes have been derived using 
Newtonian aerodynamics in conjunction with linear membrane theory, 
experimental pressure distributions obtained on other bluff shapes 
indicated that Newtonian theory does not describe the actual pressure 
distribution. However, better agreement has been obtained for many bluff 
shapes with sonic corners by means of integral relation theory. Conse- 
quently, a computational procedure involving the use of a one -strip 
integral relation te chni que for calculating pressure distributions was 
devised for use in deriving new tension shell shapes. Inasmuch as the 
nose radius and the free -stream Mach number affected the integral rela- 
tion pressure distributions, their effects on integral-relation -derived 
tension shell shapes were investigated. 

Experimental pressure distributions were obtained at a Mach number 
of 3.0 in the Langley 9 - by 6-foot thermal structures tunnel and are 
presented for typical Newtonian- and integral-relation -derived tension 



shell shapes. The experimental pressures are presented in tabular and 
graphical form and are compared with pressures obtained from the New- 
tonian and integral relation theories. Tension shell shapes were then 
derived using the experimental pressure distributions in conjunction 
with linear membrane theory and were compared with the Newtonian- and 
integral-relation-derived tension shell shapes. 
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VI. INTRODUCTION 


The possibility of exploring the atmospheres cf certain planets has 
stilulated the development of aerodynamic decelerators (refs. 1 through 
15 )• Die entry of unmanned probes into an atmosphere requires a vehicle 
with low structural weight that can generate a high drag coefficient in 
order to obtain a reasonable dwell time in the atmosphere for the gather- 
ing and transmission of data. One type of decelerator that shows promise 
of satisfying the above requirements is the tension shell as derived in 
reference 1. A typical tension shell configuration is shown in figure 1. 
The payload is assumed to be concentrated in the forward portion of the 
vehicle which is connected to the compression ring at the rear by a thin 
shell or membrane whose shape is characterized by negative Gaussian 
curvature. The configurations obtain their name from the fact that the 
shell is shaped so as to resist the aerodynamic loading by means of 
tensile stresses only. 

As indicated by the references 2 through 10, the tension shell con- 
cept has been extensively studied aerodynamically over a wide range of 
Mach number, Reynolds number, and fineness ratios. Hie results from the 
investigation of reference 2 indicate that the relatively bluff tension 
shell shapes - that is, those that generate a detached bow shock wave - 
may be best suited for use as entry vehicles because they can develop 
relatively high drag coefficients without the adverse boundary -layer 
flow separation effects that are encountered by less bluff tension shell 
shapes. In previous investigations of this concept, the tension shell 
shapes were derived using linear meribrane theory in conjunction with 
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either an axisynmetric Newtonian or a constant pressure distribution. 
However, the experimental data of reference 1 6 for large angle cones 
indicated that the constant pressure distribution predicted by Newtonian 
theory does not describe the actual distribution to which the cones are 
subjected. Thus, this inadequacy of the Newtonian theory with respect 
to cones having proportions similar to those of the bluff tension shell 
shapes poses the question of how the shapes of reference 1 might differ 
from shapes that are derived using a more realistic pressure distribution. 
Moreover, the question arises concerning the nature of the actual pres- 
sure distribution on bluff tension shell shapes inasmuch as no experi- 
mental pressures are presently available. 

In view of the above questions, an investigation was made to deter- 
mine a more realistic pressure distribution and its effect on the bluff 
tension shell shapes of reference 1. One theoretical approach that 
accurately predicts the pressure distribution of large -angle cones is 
the integral relation theory explained in reference l6. Consequently, 
a computational procedure was devised in the present investigation 
involving the use of a one -strip integral relation technique for calcu- 
lating pressure distributions and linear membrane theory to derive new 
tension shell shapes. Shapes were also derived using pressures obtained 
from wind-tunnel tests conducted at Mach 3*0 on models with shapes that 
were derived from Newtonian pressures and from pressures predicted by 
the integral relation method. 

The present thesis presents the analysis for shape determination, 
discusses the results from the theoretical study and from experiment. 
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and compares the shapes obtained from the different pressure distri- 
butions. Inasmuch as the pressure distributions given by the integral 
relation method cure somewhat dependent on Mach number and nose radius, 
the effects of varying the Mach number from 2.5 to 7.0 and of varying 
the nose radius from 0.05 to 0. 56 times the base radius on the derived 
shapes are also discussed. 



VII. ANALYSIS 


A. Statement of Problem 

The following analysis consists of solving two sets of independent 
differential equations with their respective boundary conditions. The 
first set of equations cosies from the linear membrane equilibrium equa- 
tions for a shell of revolution subjected to an axisymnetric pressure 
distribution. Circumferential and meridional stress resultants are 
assumed related to each other by a constant, and zero axial forces axe 
assumed cm the compression zing. The second set of equations is derived 
by applying integral relation theory to the solution of supersonic, 
inviscid flow around bluff bodies. The body surface is assumed to be 
normal to the free-stream flow direction at the nose, and the local 
velocity is assumed to be sonic at the sharp corner on the base com- 
pression ring. The two sets of equations are related in that the first 
'set of equations requires a pressure distribution as an input in order 
to provide a shape, while the second set of equations requires a shape 
as an input in order to provide a pressure distribution. Either a 
particular shape or a pressure distribution must be assumed to start 
the calculations, and an iterative procedure must be followed between 
the two sets of equations until a unique shape and pressure distribution 
are obtained. 


B. Basic Equations and Solutions 
1. Structural . - For a shell of revolution subject to an axisym- 
metrlc pressure distribution, the appropriate linear membrane equilibrium 
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equations are (see ref. 1): 
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Assuming N a = aliq,, where a is a constant, equation (l) becomes 


a + i 

dr 


- 0 


( 5 ) 


Solving for N^: 


= constant = NqT^ 1 " 0 '^ (6) 

Let Z (which corresponds to the value of A 2 used in ref.l) be 
given as follows: 



( 7 ) 
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With the use of equations (3), (4), (6), and (7), equation (2) becomes 


dSP + a sin q> 2P / r \^~ a ^ 

dr r cos cp cos cp \r^ J 


If R = -L, then 
*b 



and equation (8) becomes 


dflp + a sin cp _ Zp , R *(l-g) 
dR R cos <P cos <p 


(8) 


(9) 


Assuming that there is no axial force on the compression ring, the slope 
of the model surface becomes zero at r = r t , or 


dx 

dr 


-tan <p 


0 


( 10 ) 


Thus, for any particular pressure distribution , there exists a first- 
order differential equation and the necessary boundary conditions to 
solve for the model surface slope. Equation (9) may again be integrated 
with respect to r to obtain the x -coordinate of the desired configura- 
tion as a function of r. 

Equation (9), subject to the boundary condition (eq. (10)), has 
been solved for a body of revolution with an axi symmetric pressure dis- 
tribution by the use of a digital computer. A source program in 
Fortran IV language is given in appendix A. Equation (9) was numerically 
integrated twice using the fourth-order Runge-Kutta integration techni- 
que* The integration was started at R = 1 and continued with decreas- 
ing values to R = 0. In order to use the program, values of Z, a, 
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integration step height, DR, and a pressure distribution as a function 
of R are necessary inputs to the program. The program output consists 
of values of the axial coordinate and surface slope as a function of R. 
With this program, shapes axe obtained with zero nose radius of 
curvature. 

2. Aerodynamic . - For a blunt body of revolution with sharp 
corners subjected to inviscid supersonic flow, the governing differential 
equations for the one -strip method may be written as follows (see ref. 16): 

~ = (1 + KB) tan A (11) 


dp 

ds 


(1 + 


KB)^tan A - 


6 sin 6 


r 0 
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p i u i v i " ( 2 + ^'7” 
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1 V 1 


KBPqUq 
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1 + KB + 




dPl^lV 1 

* ) 


( 12 ) 


and 



A sketch of the geometry and coordinates are shown in figure 3* A 
particular body contour is specified by giving the surface angle and 
curvature as a function of s. On the axis of symmetry at s = 0, the 
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body surface must be normal to the stream direction 0(0) = £ and the 
surface slope must be continuous; thus, the following conditions hold: 

0(0)= | (Ik) 

and 

u(0) = 0 (15) 

The surface speed is required to reach sonic velocity at the model 
corner, r = r^, which results in the boundary condition that 

Uq(s*) = a* (16) 

where a* is a constant dependent on and y. Thus, we have three 

interconnected first-order differential equations (11 ), (12), and (13), 
and three boundary conditions, (l4 ), (15 ), and (l6), that must be satis- 
fied. Hie functions at the shock wave are explicit functions of y , 

3, and 0. The main dependent variables sure 5, p, and uq, while p 0 
and p Q are obtained as explicit functions of y , M^,, and Uq by using 

the isentrppic law. 

Equations (U), (12), and (13) and the boundary conditions (l4), 
(15), and (l6) have been programed for use on a digital computer. The 
integration of equations (11 ), (12), and (13) starts at s = 0 and 
terminates at s = s*; the initial shockwave standoff distance 5(0) is 
unknown and must be chosen so that equation (l6) is satisfied at the 
correct comer location. A discussion of the techniques used in solving 
the equations and a detailed program list-out are given in reference l6 
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with sample calculations on certain blunt axisymmetric bodies. In order 
to use the ex ib ting program for more general shapes, it is necessary- 
only to modify the input statements and the subroutine called BSR. In 
order to perform the calculations, an input statement was added to 
direct the computer to read in tabulated values of surface slope and 
r for any particular configuration. The subroutine BSR was modified 
by adding appropriate statements to direct the computer to use the 
correct values of surface slope and r for the particular configuration 
under consideration. With these modifications, flow conditions may be 
computed for any bluff body for which the surface slope is continuous 
and is given as a function of r. 

C. Calculation Procedure 

In order to begin the iterative procedure involving equations (9) 
and (10 ) and equations (11) through (16), either a shape or a pressure 
distribution must be assumed. For the work contained herein, an 
initial shape was calculated using a Newtonian pressure distribution 
(i.e., P = Cp = 2 sin^G ), and this shape was used in the integral 
relation computer program involving equations (ll) through (l6) to 
obtain a new pressure distribution. To satisfy the requirements of the 
integral relation computer program that the surface slopes be continuous 
and normal to the stream direction at the nose, the initial shape was 
given a spherical nose radius. The pressure distribution so obtained 
was fed into the linear membrane program involving equations (9) and 
(10) which calculated a new shape having a zero nose radius of curvature. 
The new shape was then spherically blunted and applied to repeat the 
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procedure until convergence occurred, and a unique shape and pressure 
were obtained. The process was considered to have converged if the 
difference between successive iterations resulted in a maximum variation 
in x/r^j of less than 0.0001. 



VIII. EXPERIMENT 


A. Models 

Two pressure distribution models were wind tunnel tested as part of 
this investigation. The shapes and pertinent model dimensions are given 
in figure 4; model coordinates and orifice locations are given in 
table I. The model shapes were derived using the linear membrane equi- 
librium equations for values of Z = O .65 and a = 0. For one of the 
shapes, a Newtonian pressure distribution was assumed, and for the other 
shape a pressure distribution predicted by integral relation theory was 
used. Both shapes hod a nose -radius -to-base -radius ratio (r n /r ^ ) of 
0.20. Each mode] 1 was instrumented with pressure orifices; bl orifices 
were distributed along the front face of the model and 8 orifices were 
distributed along the model base. The orifices along the frcnt face of 
the models were positioned along two radial lines l 80 ° apart and were 
mounted flush with, and normal to, the model surface, whereas the base 
pressures were measured at the open ends of tubes soldered along the 
model base (see fig. 5)- The models were machined from mild steel, and 
the surfaces were polished to a smooth bright finish. 

B. Test Apparatus 

1. Test facility . - All tests were conducted in the Langley 9“ by 
6-foot thermal structures tunnel. This facility is a supersonic blow- 
down wind tunnel which operates at a Mach number of 3-0 at stagnation 
pressures from 50 to 200 psia and at stagnation temperatures from ambient 
to 2,000° F. The air storage capacity is sufficient to permit tests of 
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2 minutes* duration for stagnation pressures of 50 psia. The models were 
sting mounted as shown in figure 6. The models were aligned at zero 
angle of attack with respect to the tunnel walls. 

2. Instrumentation . - Surface pressures were measured on each of 
the configurations by means of pressure transducers located outside the 
tunnel. The transducers were connected to the orifices by approximately 
25 feet of 0.090-inch-inside-diameter steel tubing. In addition to the 
surface pressure orifices, four total pressure and four total temperature 
probes were mounted on the walls of the tunnel to monitor the free-stream 
flow conditions. Transducers with accuracies cf ±1 percent of the maximum 
range were used for all pressure measurements . Care was exercised in 
choosing transducers having a maximum range that matched the pressures 
to be measured as closely as possible. The outputs from all the pressure 
transducers and the thermocouples were recorded by the Langley central 
data recording facility. 

C. Test Procedure 

All tests were conducted at a stagnation temperature of 250° F and 
at a stagnation pressure of 60 psia. The corresponding Reynolds number, 
based on the maximum body diameter, was approximately 10. ^ x 10^. 

Constant flow condictions were maintained for approximately 40 seconds 
in order to be certain that all pressures had stabilized. Both model 
shapes were tested twice at identical flow conditions in order to check 
for incorrect pressure readings and to determine the experimental 
accuracies involved. For the repeat test, the pressure transducers were 
interchanged. 



DC. THEORETICAL RESULTS AND DISCUSSION 


A* Comparison of Newtonian and Integral Relation 
Pressure Distributions and Corresponding 
Tension Shell Shapes 

Tension shell shape coordinates that were computed using Nevtonian 
pressure distributions and the corresponding coordinates that were 
obtained from the iteration procedure using integral relation theory- 
described previously are presented in table II. For the present work, 
three Newtonian-derived tension shell shapes were used as initial shapes 
in the iteration procedure and were computed for values of the shape 
parameter, Z, of 0 . 50 , O. 65 , and 0.80. As the shape parameter is 
increased, the body length increases. Consequently, these values of Z 
were chosen to yield shapes that were sufficiently bluff to generate a 
detached bow shock wave and thus permitted use of the integral relation 
theory. For these shapes, a nose radius of r n = 0.05 r^ was used in 
the integral relation computer program to obtain the pressure distribu- 
tions used in the linear membrane computer program for determining 
shapes. Also included in table II are the shape coordinates that 
resulted from the studies made to observe effects of nose radius and 
Mach number on tension shell shapes derived from integral relation 
theory. In both of these studies, the shape for Z = O .65 was used. 

For the study on nose radius effects, shapes were computed for nose 
radii of r n = 0.20 r^ and 0 . 5 6 r^ and for M® = 3-0. For the study 
on Mach number effects, shapes with r n = were computed for 

Mach numbers of 2 . 50 , 5*00, and 7 . 00. Corresponding pressure 
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distributions and drag coefficients predicted by the integral relation 
theory are listed in table III for all of the shapes given in table II. 

The disparity between pressure distributions given by different 
theories for the same shape is illustrated in figure 7* Presented are 
the pressure distributions calculated from Newtonian, modified Newtonian, 
and integral relation theories for the Newtonian -derived tension shell 
shape for Z = O.65 and r n = 0.05 r^. 'Hie results are typical for all 
shapes considered herein and indicate considerable difference in pre- 
dicted trends. For example, pressures obtained from the integral rela- 
tion theory decrease from the stagnation point and are substantially 
higher over most of the surface with respect to the Newtonian and modi- 
fied Newtonian values, whereas the Newtonian and modified Newtonian 
pressures increase from the nose cap tension shell juncture. The modified 
Newtonian pressures are presented in order to compare theoretical pres- 
sures based on the same free-stream Mach number; Newtonian theory is 
based on = 00 • Except for the stagnation point values, the pressures 

obtained from modified Newtonian theory (Cp = Cp,stag sin 2 ©) are in no 
better agreement with the pressures given by integral relations theory 
than are the Newtonian pressures (Cp = 2 sin 2 ©). 

The initial Newtonian -derived tension shell shape, the final 
iterated shape, and the pressure distributions used in their derivation 
are shown in figure 8 for Z = 0. 50, O. 65 , and 0.80. The final iterated 
shapes shown were obtained after four iterations. This rapid convergence 
is demonstrated in figure 8(b) for the Z = O.65 shape. Typically, a 
substantial decrease in the overa’l length of the shape occurred with 
the first iteration. Second, third, and fourth iterations resulted in 
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consecutively smaller changes in overall length as shown by the detail 
in figure 8 (b). Each successive iteration resulted in a shape with an 
overall length that was alternately shorter or longer than the preceding 
shape and thus indicated convergence. The pressure distributions con- 
verged in a similar manner, but smaller differences were obtained between 
successive iterations. The integral relation pressure values produced by 
the first and fourth iterated shapes are almost identical and are only 
slightly larger than the integral relation pressures calculated for the 
original Newtonian shape. 

The results of figure 8 show that the pressures obtained from 
integral relation theory decrease as the tension shell length is 
increased but the changes in pressure are small compared to the changes 
in shape. In contrast, the Newtonian pressures show a strong sensitivity 
to shape change. 

B. Effects of Nose Radius and Macu Number on Shapes and 
Pressure Distributions Derived From 
Integral Relation Theory 

Since the integral relation pressure distributions and the derived 
tension shell shapes are somewhat dependent on the assumed nose radius 
and the free -stream flow conditions, it is desirable to document the 
effects of these variables. Therefore, pressure distributions and 
tension shell shapes were computed for L = O .65 at a Mach number of 
3.0 with r n /rfc = 0.20 and O.56 and for Mach numbers of 2 . 50 , 3-00, 

5.00, and 7.00 with r n /r^ = 0.05* The effects resulting from the nose 
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radius and Mach number variations are shown in figures 9 and 10 , 
respectively. The largest effect of nose radius on the pressure distri- 
bution occurs over the spherical nose cap and is characterized by a 
slight bulge which covers an increasing percentage of the surface area 
as the nose radius is increased, but the overall effects are small and 
substantiate the conclusions presented in reference 17 for configurations 
with the sonic point located at the shoulder. Therefore, tension shell 
shapes that are computed from the integral relation pressure distribu- 
tions are insensitive to nose radius and a spherical nose cap may be 
added without appreciably affecting the shape downstream of the nose cap 
juncture. In contrast to these results, figure 10 shows that Mach number 
affects both the pressure distribution and the tension shell shape. An 
increase in Mach number results in a reduced static pressure loading and 
a shorter tension shell shape, but both the pressure distribution and the 
tension shell shape converge as the Mach number increases. The Mach 
number dependency and the fact that aerodynamic decelerators will be used 
over a finite Mach number range indicate that a tension shell shape 
derived for a specific application will be a compromised shape taking 
into account Mach number effects. Although a Mach number dependency is 
shown for the tension shell shape, the changes in shape that occur within 
the Mach number range considered herein are not as great as those 
obtained between Newtonian-derived and integral-relation-derived tension 
shell shapes. 



X. EXPERIMENTAL RESULTS AND DISCUSSION 


A stannary of the experimental pressures obtained from the Mach 3.0 
wind-tunnel tests is provided in table IV in pressure-coefficient form 
and is presented graphically in figure 11 for the Z = 0.65 Newtonian- 
and integral-relation-derived tension shell shapes with r n /r^, = 0.20. 

The data from both models show excellent repeatability and indicate 
nearly identical pressure distributions in spite of the difference in the 
shapes. (See fig. 4. ) Urns, the experiments! results substantiate the 
conclusion, noted earlier from the computed integral relation pressure 
distributions of figure 8, that the pressure distribution for a bluff 
tension shell shape is relatively insensitive to significant changes in 
the tension shell shapes.. For these shapes relatively high and nearly 
constant pressures are generated along the front surface to a value of 
r/r^ =» 0.75 before showing the influence of the flow expansion around 
the sharp corner at the base. Nearly constant values of pressure less 
than free- stream static pressure are obtained along the rear surface. 
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XI. COMPARISON OF THEORY AND EXPERIMENT 


A. Pressure Distribution 

In figure 11, curves of the pressures obtained from integral 
relation, Newtonian, and modified Newtonian theories are included, and 
it is apparent that the experimental data favor the curves given by- 
integral relation theory. The agreement between experiment and integral 
relation theory is excellent over the nose cap and is within 9 percent 
over the remaining portion of the Newtonian-derived shape (fig. 11(a) ) 
and within 7 percent for the iterated shape (fig. 11(b)). The 
experimental 1 y determined pressures lie above the integral relation 
curve downstream of the nose-cap juncture, and the maximum deviation 
from theory occurs at a value of r/rj, « 0.90. The deviation of integral 
relation theory from experiment may be due partially to viscous effects 
which were not taken into account by integral relation theory. These 
viscous effects result in a buildup of the boundary-layer displacement 
thickness along the model surface and therefore change the effective 
shape. However, integral relation theory gives a much better representa- 
tion of the aerodynamic loading on bluff tension shell shapes than do the 
Newtonian theories and is reccnmended for use in design applications. 

B. Shapes Determined From Experimental Pressures 
The experimental pressure distributions of figure 11 were used in 
the linear membrane computer program to determine shapes for comparison 
with the theoretically determined shapes. The coordinates of the 
experimental shapes and their corresponding experimental pressure 
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distribution in terns of the nondimensional msmbrane pressure differen- 
tial, P, are given in table V. The pressures are an average of those 
obtained from the two tests conducted on each wind-tunnel model. These 
experimental shapes and their pressure distributions a-e compared with 
corresponding Newtonian-derived and integral-relation-de.-ived shapes and 
pressure distributions in figure 12. In this figure, unlike figure 11, 
the experimental pressure values are generally less than obtained from 
integral relation theory because the experimental values of the 
parameter, P, were evaluated using measured values of base pressure, 
whereas the values of P determined from integral relation theory were 
evaluated assuming zero base pressure. The agreement between the experi- 
mental and integral relation values of P, however, is within 4 percent. 
The tension shell shapes derived from the experimental pressures differ 
by less than 2 percent and should be representative of shapes obtained 
under true aerodynamic loading. Moreover, these shapes are in excellent 
agreement with the integral- relation-derived shape. Consequently, a 
truer representation of a bluff tension shell shape can be obtained from 
pressures determined from integral relation theory rather than from 
Newtonian theory. 


C. Drag Coefficients 

One criterion governing the final selection of a planetary 
atmosphere-entry decelerator shape is the drag coefficient. Therefore, 
the method used in estimating the drag coefficient is of paramount 
importance. In table VI, the drag coefficients predicted by Newtonian- 
and integral-relation theories for the Z = 0.65 tension shell shapes 
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with r n /r^ = 0.20 are compared with drag coefficients obtained by 
integrating the experimental pressure distributions of figure 11. All 
drag coefficients are based on a free-stream static base pressure in 
order to provide a common basis for comparison. The drag coefficients 
predicted by Newtonian theory are up to 13 percent greater than experi- 
ment, whereas the drag coefficients predicted by integral relation theory 
are 6 percent less than experiment. Thus, conservative and more accurate 
estimates are obtained from integral relation theory. 



XII. CONCLUDING REMARKS 


A theoretical and an experimental investigation were conducted for 
the purpose of improving accuracy in calculating pressure distributions 
for various tension shell decelerator shapes and to evaluate the effects 
on derived shapes resulting from the application of different pressure 
distributions. Although tension shell shapes have been derived using 
Newtonian aerodynamics in conjunction with linear membrane theory, 
experimental pressure distributions obtained on other bluff shapes 
indicated that Newtonian theory does not describe the actual pressure 
distribution. However, better agreement has been obtained for many 
bluff shapes with sonic corners by means of integral relation theory. 
Consequently, a computational procedure involving the use of a one-strip 
integral relation technique for calculating pressure distributions was 
devised for use in deriving new tension shell shapes. The results 
indicated that tension shell shapes that are derived using pressure 
distributions predicted by integral relation theory are substantially 
blunter than equivalent Newtonian tension shell configurations. Although 
the pressure distributions predicted by integral relation theory are 
somewhat dependent on model nose radius and Mach number, the differences 
in the integral-relation-derived tension shell shapes attributable to 
various nose radii and Mach numbers are considerably less than those 
noted between Newtonian- and integral-relation-derived tension shell 
shapes. 

Experimental pressure distributions were obtained at a Mach number 
of 3*0 in the Langley 9- Ly 6-foot thermal structures tunnel for typical 
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Newtonian- and integral-r elation-derived tension shell shapes. The 
experimental pressure distributions were in relatively good agreement 
with theoretical pressure distributions predicted by integral relation 
theory but showed poor agreement with pressure distributions predicted 
by Newtonian and modified Newtonian theories. Tension shell shapes 
derived using the experimental pressure distributions were in good 
agreement with the corresponding integral-relation-derived tension shell 
shape but were in poor agreement with Newtonian tension shell shapes. 
Thus, integral-relation-derived tension shell shapes are more representa- 
tive of bluff tension shell shapes than Newtonian-derived shapes. Drag 
coefficients predicted by integral relation theory were approximately 
6 percent less than experimental values whereas Newtonian drag coeffi- 
cients were up to 1? percent higher than experiment. 
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XV. APPENDIX A 


In this section the Fortran IV source program is listed for the 
derivation of tension shell structures of revolution subject to linear 
membrane theory and axisymnetric pressure distribution. The following 
definitions are used in the program and are defined below or in the list 
of symbols: 


Program 

Z 

ALPHA 

THETAD 

THE(J) 

P(J) 

J 

K 

DR 

THET 

RO 

YO 


Definition 

Z 

a 

Limiting configuration slope, deg 
Radial position array, R 
Pressure dis . 'ibution array, P 
Indexing Parameter 
Indexing Parameter 
Integration step size 
0, deg 
R 
X 
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PROGRAM SAWYER < INPUT , OUTPUT ,TAPE5« I NPUT , TAPE6-OUT PUT , PUNCH » 

C PROGRAM FOR CALCULATION OF TENSION SHELL SHAPES 

100 FORMAT! 1X2H2*E16*8*2X6HALPHA«E16*8»2X7HTHETAD S E16*8) 

101 FORMAT! 8X5HTHETA8X2HRQ8X2HY012X2HP0 I 

1 FORMAT! 3E16.8) 

2 FORMAT! 2E16.8I 
20 FORMAT! 4E 16.8) 

DIMENSION THE ! 500 I »P!500) 

13 REA0!5»llZ»ALPHAtTHETAD 

READ! 5 »2 I (THE!J)?P!J) » J»l*500l 

DEG»57. 295780 

THETAS*!. 5707962 

THETAR*THET AO/DEG 

RAN*1. 5707000 

J*1 

RO*THEU> 

OR 0 s *005 
OR-ORO 

25 THETA*THETAS 
OTHET*THETA 
R=RO 
Y0*0. 

K*1 

WRITE !6, 100 )Z* ALPHA, THE TAD 
WRITE (6*101) 

THET*THET A*DEG 
WRITE!6»20)THET»RO»YO*PO 
PUNCH 2'THETAtRO 
L*1 

C COMPARISON OF ANGLES 

3 IF!R.LT*THEI JI.ANO.R.GT.THE! )G0 TO 4 

IFIR.EQ.THEC JHGO TO 5 

IFIR.EQ.THEf J+lllGO TO 6 
IF«R # LT.THE!J*1»)G0 TO 7 
IF(R.GT.THE( JH6C TO 8 
C INTERPOLATION FOR THETA A NO PRESSURE 

4 DOGMR-THE! J) l/!THE! J+ll-THEI J) ) 

PO=P(JIMP(J*l)-P(JH*DOG 

GO TO 9 

5 PO=P(J) 

GO TO 9 

6 PO=P( JA-1 1 
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GO TO 9 

7 

GO TO 3 

8 J*J-i 
GO TO 3 

C RUNG-KUTTA INTEGRATION 

9 IF(R*LE*0)G0 TO 23 
ER»R**<1. -ALPHA) 

IF (THETA. GT*RAN)A*0* 

IF (THETA. GT.RAN)C«2**2*ER*PO 
lF(THETA*GT*RAN)GO TO 14 
T AN* S IN (THETA I /COS I THE TAI 
A*-OR* ALPHA/ ( R*TAN ) 
C*2**Z*ER*PC/SIN( THETA) 

14 IF(K*EQ*2)GO TO 10 
IF(K*EQ* 3)GG TO 11 
(F(K.E0*4)G0 TO 12 
Ai*A-OR*C 

IF(THETA*GT .RAN) 91*0* 

IF ( THETA* GT • RAN) GO TO 15 
B1*DR/TAN 

15 R*R0— DR/2* 

THETA*0THET+Al/2. 

K=K*1 

GO TO 3 

10 A2*A— DR*C 

IF (THETA* GT*RAN) 82*0* 

IF ( THETA* GT* RAN) GO TO 16 
B2*0R/TAN 

16 K*K*1 

THETA»0THET*A2/2* 

GO TO 3 

11 A3=A— OR*C 

IF ( THETA* GT *RAN)83*0* 

IF (THETA *GT*RAN)GO TO 17 
B3=DR/TAN 

17 K=K*1 
R*RO— OR 

THETA=DTHET*A3 
IF (R*LE*0)C*0 
IF(R.LE*0)G0 TO 12 
GO TO 3 

12 A4=A-0R*C 
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IF(THETA.GT.RAN)94=0. 
IF<THETA.GT.RAN)GO TO 18 
B4=DR/T AN 

18 OELTA*l./6.*IAl+2.*A2+2.*A3+A4l 
OTHET =OTHET + OELTA 
THETA=DTHET 
THET=THETA*OEG 
RO=RO-OR 

OY=l • /6 • * ( B 1 *2. *62 + 2. *83+84 ) 
YO=YO+DY 

C PRINT RESULTS 

MRITE*6»20ITHET,RO»YO,PO 

PUNCH 2*THETA,RO 

K=l 

IF (R«LE* 0 ) GC TC 23 
!FtTHETA*GT.THETAR IGO TO 3 
23 CONTINUE 
STOP 
END 
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TABUS I. - COORIOVA3BS AID ORIFICE LOCATIONS FOR IEHIOIIAR- AID 

HTEGRAL-RBIATI0I-1SRIVED PRESSURE DISTRIBUTION NOEELS. 


ORIFICE 

NO. 

NEWTONIAN 

SHAPE 

■Em 

ORIFICE 

NO. 

NEWTONIAN 

SHAPE 

INTEGRAL RELATIONS 
SHAPE | 

mm 

WSM 

WSM 

Dl 

d 

d 

WSM 

BEE 

1 

0.4113 

0 

0.3588 

0 

26 

0.0811 

0.625 

0.0630 

0.625 

2 

0.4097 

0.025 

0.3572 

0.025 

27 

0.0712 

0.650 

0.0549 

0.650 

3 

0.4049 

0.050 

0.3525 

0.050 

28 

0.0619 

0.675 

0.0474 

0.675 

4 

0.3967 

0.075 

0.3442 

0.075 

29 

0.0532 

0.700 

0.0404 

0.700 

5 

0.3845 

0.100 

0.3320 

0.100 

30 

0.0450 

0.725 

0.0339 

0.725 

6 

O.368O 

G.125 

0.3167 

0.125 

31 

0.0375 

0.750 

0.0279 

0.750 

7 

0.3510 

0.150 

0.3010 

0.150 

32 

0.0306 

0.775 

0.0226 

0.775 

8 

0.3341 

0.175 

0.2854 

9.175 

33 

0.0244 

0.800 

0.0178 

0.800 

9 

0.3174 

0.200 

0.2701 

0.200 

34 

0.0188 

0.825 

0.0135 

0.825 

10 

0.3008 

0.225 

0.2550 

0.225 

35 

0.0139 

O.85O 

0.0099 

0.850 

11 

0.2845 

0.250 

0.2401 

0.250 

36 

O.OO98 

0.875 

0.0068 

0.875 

12 

0.2685 

0.275 

0.2256 

0.275 

37 

O.OO63 

0.900 

0.0043 

0.900 

13 

0.2527 

0.300 

0.2114 

O.3OO 

38 

O.OO36 

0.925 

0.0024 

0.925 

14 

0.2371 

0.325 

0.1975 

0.325 

39 

0.0016 

O.950 

0.0010 

0.950 

15 

0.2219 

0.350 

0.1839 

0.350 

40 

0.0004 

0.975 

0.0002 

0.975 

16 

0.2071 

0.375 

0.1707 

0.375 

41 

0.0001 

O.987 

0.0001 

0.987 

17 

O.I925 

0.400 

0.1579 

0.400 


0 

1.000 

0 

1.000 

18 

0.1784 

0.425 

0.1456 

0.425 

42 

* 

0.987 

* 

0.987 

19 

0.1646 

0.450 

0.1336 

0.450 

43 

* 

0.900 

* 

0.900 

20 

0.1513 

0.475 

0.1221 

0.475 

44 

* 

0.800 

* 

0.800 

21 

0.1384 

0.500 

0.1111 

O.5OO 

45 

* 

0.700 

* 

0.700 

22 

O.1259 

0.525 

0.1005 

0.525 

46 

* 

0.600 

* 

0.600 

23 

0.1140 

0.550 

0.0904 

0.550 

47 

* 

0.500 

* 

0.500 

24 

0.1025 

0.575 

0.0808 

0.575 

48 

* 

0.400 

* 

0.400 

25 

0.0915 

0.600 

0.0716 

0.600 

1 

49 

* 

0.300 

* 

0.300 


* Orifices installed on the base of the model 
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TABLE IV. - EXPERIMENTAL PRESSURE DATA FOR 'L = O .65 

TENSION SHELL SHAPES WITH r /r. = 0.20. 

a d 
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TABLE V. 




*JMEHTAL PRESSURE DISTRIBUTIONS AND THE 
DERIVED Z * O.65 TENSION SHELL SHAPES. 


^EXPERIMENTAL 
CONFIGURATION I 

^EXPERIMENTAL 
CONFIGURATION II 

P 

r/r b 

x/r b 

P 

r/r b 

x/r b 

1.830 

0 

0.3958 

1.838 

0 

0.4015 

1.815 

0.050 
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* Pressure distribution measured on Z = O.65 

Newtonian tension shell shape with r Q /r v - 0.20 

in M„ = 3.0 air stream. 

00 

** Pressure distritation measured on Z = O.65 

Integral relation tension shell shape with 

r /r, =0.20 in M = 3*0 airstream. 
n' b 00 
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TABU Ti. - THEORETICAL AID EXPHUMEITAL IRAG CQKFFICIEITS FOR 
THE Z = O .6 3 HEWTOKAI- A>v> UTEGRAIr-RElATIOB- 
BBtlVKD TEBSIOK m SHAF& WITH r /r fe * 0.20. 
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Figure 2.- Tension shell shape and coordinate system. 



Figure - Geometry and coordinate system for aerodynamic 

consideration. 














Integral Relation Theory 
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Fig-are 7.- Theoretical pressure distributions for the Z = O.65 
Newtonian-derived tension shell shape with r n /r-^ = O.05. 



Figure 8.- Comparison of Newtonian and integral relation pressure 
distributions and corresponding tension shell shapes. 
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Figure 9*“ Effects of nose radius on the integral relation pressure 
distributions and derived tension shell shape*, at H» = 3*0. 






Pressure Distribution* 





Figure 10.- Effects of Mach number on the integral relation pressure 
distributions and derived tension shell shapes with r^r-^ = 0.05 



(a) Newtonian shape. 


Figure 11 .- Experimental and theoretical pressure distributions about 
the Z = O.65 Newtonian- and integral-relation-derived tension 
shell shapes with r n /*b * 0,20 8114 at = 5*0* 
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(b) Integral relation shape 
Figure 11.- Concluded. 








